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Abstract: We present a pedagogical review of the weak gravitational lensing measurement pro¬ 
cess and its connection to major scientific questions such as dark matter and dark energy. Then 
we describe common ways of parametrizing systematic errors and understanding how they affect 
weak lensing measurements. Finally, we discuss several instrumental systematics and how they 
fit into this context, and conclude with some future perspective on how progress can be made in 
understanding the impact of instrumental systematics on weak lensing measurements. 
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1. Motivation 

In our current cosmological model, only 4% of the Universe consists of the baryonic matter from 
which stars and planets are made. To explain a wide variety of cosmological observations, we have 
been forced to posit the existence of dark matter (detected through its gravitational attraction) and 
dark energy (which causes a repulsion that is driving the accelerated expansion of the Universe). 
While we infer the existence of these dark components, the question of what they actually are 
remains a mystery. 

Gravitational lensing flU 0, ED is the deflection of light from distant objects by the matter along 
its path to us - all of the matter, including dark matter. Lensing measurements are thus directly 
sensitive to dark matter. They also allow us to infer the properties of dark energy 0], because 
the accelerated expansion of the Universe that it causes directly opposes the effects of gravity, 
which tends to cause matter to clump into ever larger structures. This measurement relies on the 
small but spatially coherent distortions (known as weak shears) in the shapes of distant galaxies, 
which provide a statistical map of cosmological large-scale structure along the line-of-sight. A key 
problem for weak lensing measurement is that there are many other effects in astronomical images 
that cause coherent distortions of galaxy shapes and that must be effectively removed in order to 
reliably measure weak lensing. 

Because of its sensitivity to dark matter and dark energy, ever-larger surveys have been planned 
to measure weak lensing increasingly precisely, including Euclid 1 [0, LSST 2 [0, and WFIRST- 
AFTA 3 [0], all of which are Stage IV dark energy experiments according to the Dark Energy Task 

J http://sci.esa.int/euclid/, http://www.euclid-ec.org 

2 http://www.lsst.org/lsst/ 

3 http://wfirst.gsfc.nasa.gov 
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Force [0] definitions. The very high signal-to-noise ratio of the weak tensing measurements in 
these surveys will necessitate excellent understanding and control of systematic uncertainties. 

Broadly speaking, among the sources of systematic uncertainty in weak lensing measurements 
are (1) instrumental systematics, (2) the need to remove the point-spread function (PSF) which 
blurs the galaxy images to robustly infer the lensing shear, (3) the need to estimate redshifts for 
the galaxies or at least statistically understand their distribution, and (4) theoretical uncertainties 
such as the impact of baryons on the measured quantities (usually predicted using dark-matter-only 
simulations) and failures in the assumption that galaxy shapes are randomly oriented in the absence 
of lensing. Of these, issues (2) through (4) have been thoroughly explored in the literature (see, for 
example, [@, HO, H3L 13, Q Q3L O, HO, O] and references therein). Instrumental systematics have 
received somewhat less attention, and are the main topic of discussion in this pedagogical review. 

We begin in section 0 by discussing how weak lensing measurements are made, and how those 
measurements are affected by systematic errors in general. In section 0 we discuss instrumental 
systematics more specifically. We conclude in section f| with future outlook. 

2. Scientific background 
2.1 Weak lensing measurements 

In this subsection, we ignore the question of systematics, and give a basic summary of how weak 
gravitational lensing measurements are carried out. 

Cosmic shear, the measurement of weak lensing by large-scale structure, relies on correlating 
pairs of galaxy shapes as a function of their separation on the sky. Galaxies that are nearby on the 
sky tend to have correlated shapes due to having been lensed by associated structures that subject 
them to similar weak lensing shears. Galaxies that are farther apart on the sky have much less 
correlated weak lensing shears. 

To carry out this measurement, the most basic requirement is a catalog of galaxy positions and 
their shape estimates 4 , which we take as a proxy for shear (but see section ^ 2 | for more details on 
the difference between galaxy shapes and shear). While galaxy light profiles do not in general have 
elliptical isophotes, they are nonetheless commonly described in terms of some effective ellipticity 
magnitude |el, which can be defined in terms of the minor-to-major axis ratio, and position angle 
0 with respect to some coordinate system on the sky. Alternatively, one can write two components 
of the ellipticity as £\ = |e|cos 20,£2 = |e|sin 20 , and define a complex ellipticity as £ = £\ + 
i£ 2 - These trigonometric functions of 20 are very common in mathematical expressions related to 
galaxy shapes or lensing shears, due to their being spin -2 quantities (ellipses are invariant under 
180° rotations). For more details on common definitions of galaxy shapes, see for example lfT911 . 

The galaxy shapes will have some large random intrinsic components which is a key source 
of statistical error in most weak lensing measurements, along with a small but fortunately coherent 
weak lensing shear. The intrinsic shape then averages out of correlation functions of galaxy shapes, 

technically, one need not have a per-galaxy shape or shear estimate to reconstruct the lensing shear field; see IfTRII for 
options that avoid calculating a per-object shear in favor of ensemble estimates. However, these methods are relatively 
new and their use for cosmic shear requires further practical and theoretical development, so for the remainder of this 
review we assume that weak lensing measurements will use per-object galaxy shapes. 
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while the tensing shear does not (but see lfT711 for issues related to this basic assumption). An esti¬ 
mate of the galaxy redshift distribution is also needed in order to make a theoretical prediction for 
the cosmic shear correlations. However, even more promising is the so-called cosmic shear tomog¬ 
raphy, which involves dividing the galaxies into redshift bins based on a per-galaxy photometric 
redshift estimate (see, for example, figure 1 in 1 1200 . With > 1 redshift bin, additional information 
can be obtained by correlating galaxies both within and across those redshift slices, and seeing how 
the shear signal changes in each case. 

A typical cosmic shear correlation function estimator, for the correlation function correspond¬ 
ing to redshift slices i and j is 11201 




LWgWb [ei(xa)ej{x b ) ± 4 (x q )£' < (Xfr)] 

'LWaWb 


( 2 . 1 ) 


The index a runs over all galaxies in redshift slice i, while the index b runs over all galaxies 
in redshift slice j. Galaxy pairs (a : b) are identified, and their angular separation on the sky 0 
(calculated using their positions x a ^) is used to put them into bins in 0. Each galaxy has a weight 
factor w a ^ that relates to uncertainty in the shear estimates. Each galaxy’s shear estimate £, which 
has two components, is rotated into a coordinate system defined by the vector connecting galaxies 
a and b. The component corresponding to orientation along or at 90° with respect to that vector 
is £*, and the component corresponding to orientation along the ±45° rotated system is £ x . The 
two estimated correlation functions for each pair of redshift bins i and j, £±(0), both relate to 
integrals over the matter power spectrum along the line-of-sight. For examples of recent cosmic 
shear measurements using several different datasets, see EDL EE EE E3h These have been used to 
constrain the amplitude of the matter power spectrum and its evolution with time, which tells us 
about the equation of state of dark energy. 

Another type of weak lensing measurement is cluster-galaxy or galaxy-galaxy lensing, which 
involves identifying specific foreground galaxies or clusters that act as “lenses”, and correlating 
them with the positions of background sources. This measurement is therefore a (lens) position vs. 
(source) shape cross-correlation function, which makes it less subject to certain systematic errors 
in shear estimates. It provides an estimate of the projected total matter surface density around 
the chosen lenses. While it has primarily been used to make the connection between galaxies 
or clusters and their host dark matter halos (e.g., very recently, EH E3D, it can also be used in 
combination with estimates of galaxy clustering to constrain the growth of cosmic structure and 
therefore dark energy EHh and to test whether General Relativity is the effective theory of gravity 
on cosmological scales EE EBD- 

Examples of both types of lensing observables are shown in figure |T|. 


2.2 The impact of systematics in shear estimation 

We now consider the main sources of systematics in estimating the lensing shears, how they are 
commonly described, and how they affect the measured statistics used for many cosmic shear 
studies in equation EH] . 

First, it is important to bear in mind that the concept of a galaxy shape is not very well- 
defined. Even moderate-resolution imaging reveals that galaxies do not have elliptical isophotes, 
but rather may have some degree of dumpiness in their light profiles and/or ellipticity gradients 
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SDSS DR7 LRGs, 0.16<z<0.36 



Figure 1 . Left: A galaxy-galaxy lensing measurement from the Sloan Digital Sky Survey (SDSS, EU) il¬ 
lustrating the projected mass profile around Luminous Red Galaxies (LRGs), including the host dark matter 
halo below l/z _1 Mpc and contributions from large-scale structure. At the largest scales, the signal corre¬ 
sponds to shears (shape distortions) of ^ 5 x 10~ 4 . Right: A theoretical cosmic shear correlation function 
(figure 14.4 from 0]), illustrating the typical magnitude and scaling of the shear correlation function with 
separation, as well as an expected signal due to one particular systematic error. 


(e.g. a rounder bulge component in the inner region and a more extended and flattened disk). Thus, 
different methods of estimating galaxy shapes will often give quite different answers for the same 
galaxies. Even the same method may give a different answer for a single galaxy when measuring 
in different passbands, due to different components of the galaxy having different spectral energy 
distributions. These differences are, from a weak lensing perspective, not of interest. What is 
important is only that a given shape estimator have a well-defined response on average to lensing 
shears, such that we can use the ensemble averages to measure statistics of the lensing shear field. 

Images from telescopes do not directly show us the shapes of galaxies. The galaxy light 
profiles are modified by many effects. First and foremost is the point-spread function (PSF) of 
the atmosphere (for ground-based telescopes) and the telescope optics. This can be modeled as 
a convolution kernel that rounds galaxy shapes by a significant factor, but can have some small 
ellipticity that gets coherently imprinted into galaxy shapes and mimics a cosmic shear signal. 
There are also distortions of coordinate systems that can be treated as affine transformations rather 
than convolutions. Images contain defects such as cosmic rays, bleed trails, and others. Finally, the 
images have noise, which results in significant biases in estimates of shear E2L EDD- A schematic 
illustrating the image generation process is shown in figure Note that the process of removing 
the PSF to infer the shears robustly has been the subject of focused attention from the weak lensing 
community for nearly two decades (with several large community-wide blind challenges [0, ITTL 
0,111 II). Nonetheless, some issues remain to be solved, even when ignoring instrumental 
systematics. 
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Figure 2. An illustration of the process of gravitational lensing and other effects that change the apparent 
shapes of galaxies in the astronomical imaging process (from EH])- Instrumental systematics are not shown, 
but would appear after the PSF convolution. 

A common description for systematic errors in shear estimates is a linear model E2D, 

7 =(l+m)y+c, (2.2) 

which can be written separately for each of the two components of the shear. Here 7 is the true 
lensing shear, 7 is the estimated shear, m is the multiplicative bias and c the additive bias in the 
estimated shears. This additive bias is in general linearly proportional to the PSF ellipticity Ifl5ll . 
though selection effects that correlate with the PSF direction might cause a deviation from this 
typical relationship. Ideally, m = c = 0. 

We can model galaxy-galaxy or cluster-galaxy lensing as a galaxy position g vs. estimated 
shear 7 cross-correlation, or (schematically) 

(gf) = (1 +m) {gy) + (gc). (2.3) 

The first term on the right-hand side is the ideal signal with a multiplicative calibration bias de¬ 
termined by m. The second, additive term relates to the correlation between galaxy positions and 
additive systematics, which is typically either zero or removable through simple means I Bill . Using 
similar notation, cosmic shear can be written as 

(77) = (1 + ra) 2 (yy) + 2(1 + m){cy) + (cc). (2.4) 

The first term on the right-hand side is the real cosmic shear with a calibration factor (1 +m ) 2 « 
1 + 2 m (since typically m « 1). The second term would involve a correlation between additive 
systematics and the true shear, which may be nonzero by chance for small regions of a survey 
but naturally averages to zero when considering large areas. Finally, the last term is a coherent 
additive term in the shear field, which can typically be diagnosed using correlations between the 
galaxy shear estimates and star shapes IK411 . Generally the scaling of this term with separation will 
depend on the correlation function of the PSF ellipticity or other source of additive systematics; for 
examples, see figure |3j. Since there are clear diagnostics for additive systematics, the weak lensing 
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Figure 3. Lensing shear correlation functions £±, along with the analogous shape correlations for PSF 
anisotropies from several sources: typical atmospheric PSFs, and space- and ground-like optical PSFs. As 
shown these can have quite different scalings with separation, and at some scales they have amplitudes above 
the cosmic shear signal, highlighting the need for very accurate removal. From IFTTII . 


community has primarily focused on multiplicative biases in the recent past due to the difficulty in 
getting an absolute calibration. 

In general, this linear model for systematic errors can be used to describe basic errors in shear 
estimation due to the process of correcting galaxy shapes for the effects of the PSF. For other 
systematics, including instrumental ones, its applicability may be limited. For example, if the 
PSF is incorrectly modeled in a way that results in the PSF being systematically too large by the 
same amount at all places on the sky, then this could be modeled as a spatially constant (m) over 
the whole galaxy population. Real errors in the modeling of the PSF tend to have some scale- 
dependence, and one has to consider the full spatial correlation function of those PSF modeling 
errors in order to really understand the impact on science. The linear model is still useful for 
providing basic intuition, however. 

It is important to note that in general, m = m(S/N, size, morphology) and therefore it is also 
effectively a function of redshift. The same is true for additive errors (c). Thus, in a tomographic 
measurement, each redshift slice must have its own effective systematic errors estimated and ac¬ 
counted for, rather than treating the entire galaxy population in a monolithic way. 


3. Instrumental systematics 

Instrumental systematics include a variety of effects such as nonlinearity, charge-transfer ineffi¬ 
ciency (most notably in the Hubble Space Telescope Advanced Camera for Surveys, (ED), tree 
rings and edge distortions fEHl . the brighter-fatter effect 115711 . image defects leading to coherent se¬ 
lection effects (masking bias, EEDX amplifier overshoot, crosstalk, fringes, and gate diffraction. At 
a lower level, there are also possible chromatic instrumental effects, which would affect galaxy and 
star images differently due to their different spectral energy distributions. Several of these effects 
(tree rings, edge distortions) can be treated as part of the WCS (world coordinate system) mapping 
from image to world coordinates, albeit a rather complex part of the WCS. They do not necessarily 
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fit cleanly into the model for systematics presented in equation [2T21 . Here we discuss two of these 
issues in more detail. 


3.1 Brighter-fatter effect 


The brighter-fatter effect ICTTI . ES is a result of charge building up in pixels and repelling additional 
electrons (pushing them toward neighboring pixels). This makes the PSF larger for brighter objects 
and smaller for fainter objects. The problem for weak lensing is that bright stars (signal-to-noise 
ratios of > 50) are typically used to estimate the PSF, while the galaxies used to estimate weak 
lensing shear are on average quite faint. Thus, their effective PSF will be smaller than the effective 
PSF model inferred from stars, and the PSF correction routine will systematically over-correct for 
the blurring effects of the PSF, resulting in overestimates of shear. 

We can use a simple model for how the PSF affects galaxy shear measurements to estimate 
the impact of the brighter-fatter effect (if left uncorrected). Let us assume that for some particular 
CCD, the linear PSF size (e.g., its FWHM) will be overestimated by 1% due to the brighter-fatter 
effect (with 2 % being the maximum effect seen in IP7I1 using a wide range of fluxes). IflOl showed 
that in the Gaussian approximation, the systematic bias in the inferred shear (m = Sf/f) due to a 
misestimate of the PSF size is 


8f 

7 



0 


ST 

T(p) 


(3.1) 


Here, is the trace of the second moment matrix of the PSF, and ST^ /T^ is its fractional 
error. If the linear size of the PSF is misestimated by 1% then ST^/T^ « 0.02. R 2 is the galaxy 
“resolution factor”, which is 0 for unresolved objects and approaches 1 for highly resolved galaxies. 
A galaxy near the resolution limit might have R 2 ~ 0.25, which results in a shear calibration bias 
m = 0.06. Moderately resolved galaxies (R 2 ~ 0.5) would have m = 0.02. For context, large future 
surveys require that the calibration bias be reduced to the level of m < (1 — 3) x 10 -3 , or at least be 
known and reliably corrected at that level [0, ED, in order for the calibration bias to not contaminate 
the dark energy constraints at a level comparable to the statistical errors. 

If the brighter-fatter effect has some directionality, then the PSF ellipticity for typical galaxies 
might likewise be misestimated. For typical methods of shear estimation, a few percent of the PSF 
ellipticity leaks coherently as an additive c term in the galaxy shears Ifl5ll . If the PSF ellipticity 
error due to the brighter-fatter effect is 0.01 (coherently), then a typical shear estimation method 
will end up with a c term of ~ 3 x 10 -4 , again slightly above the requirements for future surveys. 

Clearly both the multiplicative and additive biases that may be induced by the brighter-fatter 
effect are larger than the requirements for upcoming surveys, so this effect must be well-modeled 
and removed. 


3.2 Nonlinearity 

An example of a long-known instrumental effect is CCD nonlinearity, which can be important at 
the high flux levels observed in the bright star images commonly used to estimate the PSF (but not 
for most galaxy images). CCD nonlinearity is commonly measured in the lab and then removed, 
making it an essentially ignored effect for most weak lensing shear studies. However, if there is 
some error in the nonlinearity estimation, or if it is a poorly-tracked function of time such that 
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the original corrections cannot be used later on, this can cause errors in the “corrected” images of 
bright stars used to estimate the PSF, meaning that the PSF model used to correct the galaxy shapes 
is wrong. 

An example of this error in practice has been seen in the SDSS. As shown in lF%ll . there is 
an error in the PSF model and shear estimates in one of the six r-band CCDs, for which the most 
likely explanation is an error in the nonlinearity corrections 5 . The impact of this effect can be seen 
in figure 8 of that work, which shows a coherent systematic error in the galaxy shapes at the level 
of 5 x 10 3 . 1 12511 showed that the majority of the additive systematic errors in the lensing shear in 
the entire SDSS area are eliminated after removing the 1/6 of the data from that one CCD. While it 
seems unlikely that the nonlinearity for CCDs for future surveys will be misestimated in advance, 
if there is any change in that nonlinearity in the years needed for these surveys to be completed, 
then systematic errors of this type are possible. 

4. Conclusions 

Instrumental systematics can leak into weak gravitational lensing measurements in surprising ways. 
We have described the quantities that are measured as part of a weak lensing measurement, and gave 
examples for how two particular instrumental systematics can affect those measurements. In future, 
more work will be needed to fully characterize the various instrumental systematics that will affect 
surveys such as the LSST. There are, fortunately, several means for doing so. Laboratory tests can 
be used to understand the magnitude of the effects, which can then be included in survey image 
simulation software such as that for LSST (PhoSim: 02L S3)- Several instrumental effects have 
also been included in GalSim 6 a, an open-source image simulation software package, which can 
be used for more focused studies of particular instrumental effects without including all possible 
systematics. In either case, these simulation packages can be used to simulate images and therefore 
estimate the effects of these systematics on the correlation functions used to measure cosmic shear 
and cluster-galaxy or galaxy-galaxy lensing, and thereby estimate how the systematic errors affect 
inferences about dark matter and dark energy. The simulations can also be used (perhaps even more 
importantly) for tests of mitigation schemes. 
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